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ABSTRACT 

Tracing the cosmic evolution of the Baryonic Acoustic Oscillation (BAO) scale with galaxy two 
point correlation functions is currently the most promising approach to detect dark energy at early 
times. A number of ongoing and future experiments will measure the BAO peak with unprecedented 
accuracy. We show based on a set of N-Body simulations that the matter distribution is anisotropic 
out to ~ 150 h~ 1 Mpc, far beyond the BAO scale of ~ 100 /i _1 Mpc, and discuss implications for the 
measurement of the BAO. To that purpose we use alignment correlation functions, i.e., cross corre- 
lation functions between high density peaks and the overall matter distribution measured along the 
orientation of the peaks and perpendicular to it. The correlation function measured along (perpen- 
dicular to) the orientation of high density peaks is enhanced (reduced) by a factor of 2 compared to 
the conventional correlation function and the location of the BAO peak shifts towards smaller (larger) 
scales if measured along (perpendicular to) the orientation of the high density peaks. Similar effects 
are expected to shape observed galaxy correlation functions at BAO scales. 

Subject headings: large-scale structure of universe — methods: numerical 
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1. INTRODUCTION 

The Baryonic Acoustic Oscillations (BAO) constitute 
a characteristic feature within the large scale struc- 
ture of the Universe and can serve as standard ruler 
for constraining the properties of dark energy (e.g., 
Blake fc Glazebrookl200"llLinderll2003l:lSeo fe Eisensteir 
McDonald & EiscnstemT 12007 



2003T IWand 1200ft _ 
Seo fe Eisensteinll2007tlSeo et al.l 2008. 2009; Kazin et al.l 



2012TT 



The BAO in the baryon-photon fluid of the pre- 
recombination era imprint the sound horizon distance 
at decoupling as a typical scale in th e matter cor- 
or 
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) . These oscillations 



were detected in the cosmic microwave background 
(e.g., iPage et all 120031) and in the spatial distrib ution 
of galaxies (jEisenstein et al.|[2005t iCole et all 120051) and 
have been confirmed by a number of subsequent studies 
e.g.. iPercival et a^l2007l: ICabre fc Gaztanagal IL 'Ut)'): 



Sanchez et al.l 1200a iPercival et al.l I2010L iReid et alJ 



2010; iKazin et al.l^dloT T 

The next generation of large galaxy surveys, like the 
Panoramic Survey Telescope fc Rap id Response Sys- 
tem fPan-STAR RS.lKaiser et al.ll2002l). the Dark Energy 
Surve y (DES, iThe Dark Energy Survey Collaboration! 
120051 ). the Baryonic Oscilla tion Spectro scopic Survey 
(TiOSS. ISchlegel et al.1 l2009h . BigBOSS (ISchlegel et alJ 
120 111) , the Hobby Eberly Telescop e Dark Energy Exper- 
iment fHETDE X. iHill et al.1 120041) and the space based 
Euclid mission (jCimatti et al.l 120091 ) . will cover volumes 
much larger than current datasets, allowing for much 
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more accurate determinations of the BAO. 

In the two-point correlation function the BAO are 
visible as a uniq ue broad and quasi Gaussian peak 
(jMatsubaral 120041 ) . However, the determination of the 
shape and loc ation of the peak may be a ffected by sam- 
ple v a riance (ICabre fc Gaztahaga 2009; Martine z et al.1 
[20091: IKazin et all 120101: ICabre fc Gaztanagal I2011D . 
non linear e ffects and the bia s of the tracer galaxy 
popu l ation dSmith et al.1 120031: iCrocce fc Scocci 
2006: lAngulo et al.l 12001 ICrocce fc Scoccimarrol 1200 
' 20081). These difficulties lead 
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([201 1[) to suggest to use the zero-crossing of the two 
point correlation located at ~ 130 /i _1 Mpc as standard 
ruler instead of the peak location. Yet, several obser- 
vations do not show the theoretically predicted zero - 
crossing at all ([Martinez et al.1 120091: IKazin et al.ll2010D . 
At this stage it is unclear whether this discrepancy 
is caused by systematic effects or cosmic variance or 
whether it repr esents a challenge for the concordance 
ACD M model (iSylos Labini et all 120091 : iSanchez et all 
120091 IKazin et al.ll2010D .~ 

One basic assumption for interpretation of the BAO 
measurements is that the matter distribution is isotropic 
at the relevant scales (~ 100/t~ 1 Mpc) . In this work 
we use alignment correla tion functions ([Paz et al.1 120081 : 
iFaltenbacher et aT1l2009f ) to show that the amplitudes of 
the two-point correlation function measured along the 
orientations of the high density peaks are larger than 
those derived from spherically averaged (conventional) 
clustering analysis out to scales of > 150 /i -1 Mpc and 
discuss possible effects on measurements of the BAO. 

2. METHODOLOGY 

We use a set of large scale dark matter simulations 
which follow the evolution of the cosmic density field in- 
cluding its collapse into high density peaks. Based on 
their orientations we compute the alignment cross corre- 
lation function with the over all mass distribution out to 
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Fig. 1. — Upper panels: 3 dimensional two point cross correlation functions between FoF groups of the indicated mass and the overall 
matter density field represented by 10% random subset of all simulation particles. Solid lines show the conventional correlation functions. 
Dotted and dashed lines show the alignment correlation functions along and perpendicular to the orientations of the FoF groups, respectively. 
Lower panels: Difference between conventional and alignment correlation functions above. The shaded areas in the upper and lower panels 
give the cosmic variance. Poisson errors are negligible and are not shown here. 

~ 200 ft, _1 Mpc. In this section we briefly discuss the sim- 
ulations, the computation of density peak orientations 
and the definition of the alignment correlation function. 



2.1. Simulations 

We used GADGET-2 (ISpringell l2005h to carry out 
50 1 /i~ 1 Gpc-simulations based on the same concordant 
ACDM cosmology but different realizations of the initial 
density field. The total volume of the ensemble is almost 
70 times larger tha n the volume used for th e detection of 
the BAO signal by Eiscnst ein et al.l (|2005l ). The cosmo- 
logical parameters co rrespond to those us ed for the Mil- 
lennium Simulation (jSpringel et al.l 12005): matter den- 
sity parameter, J7m = 0.25; cosmological constant term, 
Oa = 0.75; power spectrum normalization, as = 0.9; 
spectral slope, n s = 1; and Hubble parameter, h — 0.73. 
For each run 320 3 particles are used resulting in a mass 
resolution of 2.11 x 10 12 ft _1 M Q for the dark matter par- 
ticles. A Plummer-equivalent softening of 80 /i _1 kpc was 
employed. An analysis based on a much lower force res- 
olution showed very similar results. 

2.2. Orientations of high density peaks 

Density p eaks are identifie d by using a friends-of- 
friends (FoF. lDavis eTaIT l985) group-finder with linking 
parameter b = 0.2 times the mean particle separation. In 
each realization there are about 30,000 FoF-groups with 
masses greater than 10 14 /i _1 M Q . The orientations are 
derived by evaluating the eigenvalues of the second mo- 
ment of the mass tensor. In the following the direction of 
principal axis associated with the largest eigenvalue is re- 
ferred to as the (3D) orientation of the density peak. Ob- 
servationally 3D orientations are difficult to determine, 
therefore, we also use projected orientations. These are 
determined by collapsing the mass distribution of the 
density peak along the line of sight (here equivalent to the 
z-axis) and, consecutively, computing the eigenvectors of 
inertia tensor for the two dimensional mass distribution. 
Despite the poor force resolution the orientations of the 



density peaks are expected to trace the large scale tidal 
field well (see below). 

2.3. Alignment correlation functions 

Two-point correlation functions, £(r), are a primary 
tool for quantifying th e clustering p roperties of cosmo- 
logical structures (e.g. Peebles 1980) .We employ an ex- 
tension of the two-point cross-correlation function which 
can be used to quantify spatial alignment of objects 
in one sample with respect to orientat ions assigned to 
the objects in the othe r sample (see. iPaz et al.l 120081 : 
iFaltenbacher et al.l 120091: iSchneider et al.l 1201 ID . Given 
a sample of objects in question (Sample Q), a sample 
of reference objects (Sample R) and a random sample 
(Sample TV) we can define the alignment two point cross 
correlation function 



N n QR(e,r) 
N R QK(9, r) 
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where r indicates the distance, 9, is the angle between 
the orientation of the objects in Sample Q and the line 
connecting to objects in either of the two other samples. 
Nr and Nfi are the number of points contained in the 
reference and random sample, QR{9, r) and QTZ{9, r) are 
the counts of cross pairs between the indicated samples 
for given 9 and r. The conventional correlation function 
is the average of the alignment correlation function over 
the full range of 9 values. 

Taking symmetries into account the value of 9 ranges 
from zero (along the orientation of the main object) 
to 90 degrees (perpendicular to its orientation). Thus, 
higher amplitudes of the alignment correlation functions 
at small values of 9 indicate that the reference objects are 
more likely to be aligned along orientation of the main 
objects. In contrast, higher amplitudes at larger angles 
indicate that the reference objects are more likely to be 
located perpendicular to the orientation of the main ob- 
jects. 

We employ the alignment correlation functions, more 
exactly, the average of the alignment correlation function 
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Fig. 2. — Left panel: Same as middle panel in Fig. [T] for a redshift 
projected orientations and redshift space separations (at z = 0). 

within given angular ranges, to quantify the large-scale 
alignment between the matter distribution in the Uni- 
verse and the orientation of high density peaks (Mf f > 
10 14 /i _1 M Q ). The reference samples are generated from 
10% random subsets of the overall particle distributions. 
If real-space correlations are considered we use the three 
dimensional orientation and separation vector to deter- 
mine 6. Alternatively, if observational consequences are 
discussed we use projected orientations but still 3D dis- 
tances since the effects of reshift-space distortions are 
small. 

3. RESULTS 

In this section we present our findings for the alignment 
correlation functions in real- and redshift-space and for 
redshifts z=0 and 0.6. 

Fig. [T] shows the conventional and alignment two point 
cross correlation functions (upper panels) and the resid- 
uals of the alignment correlation functions about the 
conventional correlation function (lower panels) employ- 
ing different lower mass cuts (1 x 10 14 , 5 x 10 14 and 
1 x 10 15 h~ 1 M@) for the high density peaks in the main 
sample. Solid lines show the conventional correlation 
functions. Dotted and dashed lines represent the cor- 
relation functions along and perpendicular to the orien- 
tations of the density peaks. The gray regions represent 
the cosmic variance. Interestingly, the cosmic variance 
in the lower panels is slightly smaller, this is because the 
differences between conventional and alignment correla- 
tion functions are independent of their absolute values. 
To check for systematic or numerical errors we have re- 
peated the above analysis with randomly interchanged 
orientations within each realization. In this case we do 
not detect any significant difference between the conven- 
tional and alignment correlation functions. 

The left panel of Fig. [5] displays conventional and align- 
ment correlation functions for density peaks with masses 
larger than 5 x 10 14 h~ 1 M <3 at z = 0.6. The right panel of 
Fig. [2] represents the correlation functions based on pro- 
jected orientations. Here, the angles used to select the 
pairs along and perpendicular to the orientations of the 
density peaks are computed between the projected ori- 
entations and the projected distance vectors. Distances 
are computed in redshift space using the distant observer 



of z = 0.6. Right panel: Same as middle panel in Fig. [T] but using 

approximation. In this case the differences in the ampli- 
tudes between the conventional and the alignment cor- 
relation functions are somewhat reduced but still signif- 
icant beyond the BAO scale. 

In conclusion, the alignment correlation function 
between dark matter density peaks with masses > 
10 14 /i _1 Mq and the overall matter distribution reveals 
anisotropies at scales larger than 150 /i _1 Mpc far beyond 
the BAO scale. The difference between conventional 
correlation functions and correlation functions measured 
along and perpendicular to the orientations of the peaks 
is most pronounced for the most massive peaks in real- 
space at z = 0. But it remains visible at earlier time and 
also if projected orientations are used. 

4. DISCUSSION 

The findings presented here are based on N-Body sim- 
ulations however qualitative similar effects are expected 
to shape the observed galaxy correlation functions 
at scales up to 150 ft, _1 Mpc. Potential observational 
implications are: 

• Large-scale anisotropy and structure of two point cor- 
relation functions: Our results indicate that matter is 
distributed anisotropically out to separations far larger 
then the BAO scale. At any given pair separation the 
conventional correlation function is the average of the 
alignment correlation function over the whole range of 9. 
We find that the alignment correlation function measured 
along high density peaks does not fall below zero (for 
separations < 150 /i _1 Mpc) which is counterbalanced by 
negative amplitudes at much smaller separations for mea- 
surements perpendicular to it. At the BAO scale galaxy 
two point correlation functions can be interpreted as the 
average of a more highly clustered component along the 
direction of high density peaks and a less clustered com- 
ponent perpendicular to it. 

• Shape and location of the BAO peak: At BAO scales 
the amplitudes of the correlation functions between high 
density peaks and the overall matter distribution are sig- 
nificantly higher if measured along the orientation of the 
peaks. In this case the BAO peak is composed of a hump 
on top of the declining but still positive correlation func- 
tion. The signal perpendicular to the orientations is dom- 
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inated by the B AO hump itself with negligible underlying 
clustering signal. The location of the BAO peak is found 
at somewhat smaller (larger) scales if measured parallel 
(perpendicular) to the orientation of high density peaks. 
For the measurements along the density peak orienta- 
tions the BAO hump is transformed into a plateau-like 
feature. If the survey volume is dominated by large fil- 
amentary structures the shape of the conventional cor- 
relation function should be close to that found for the 
pa rallel signal s hown here. This may relate to a study 
bv lKazin et al.l (|2010l ) who find no apparent peak in 6% 
of their mock samples. 

• Zero crossing and large scale power: Several publi- 
cations report the non-detection of zero-crossing out to 
scales of 250 fe~ 1 Mpc which indicates ( unexpected) large 
scale power (e.g-.lMartinez et al.ll2009t iKazin et alJl2010b 
iSvlos Labini et al.ll2009[ )~ We find similar results for the 
correlation functions measured along the orientations of 
high density peaks. Furthermore, the alignment cross 
correlation function measured along the orientations of 
high density peaks at z = 0.6 shows zero-crossing in con- 
trast to the behavior at z = 0. The WiggleZ redshift- 
space correlation function at z = 0.6 (jBlake et al.l 120 111 ) 
shows a crossover as well. Whether these analogies are 
coincidental remains to be explored in future work. 

• Direct measurement of anisotropy: With the ad- 
vent of enormous cl uster catalogs (e.g. JHao et al.|[2~010t 
IGilbank et al.ll20TT[) it should in principle be possible to 



directly measure the large scale anisotropies with align- 
ment correlation functions if cluster orientations can be 
determined with sufficient accuracy. An anisotropy sig- 
nal may even be extracted simply by using orientations 
of the cluster central luminous red galaxies because the 
orientati ons of central galaxies and host systems are cor- 
related (iFaltenbacher et al.l 120091: lOkumura et~all 120091: 
iSchneider et al.ll2011l) . ~ 

• Improvement of BA O measurements: If the directional 
effects reported here are observable it would be worth- 
while considering to measure the BAO peak perpendicu- 
lar to the orientations of galaxy clusters (or luminous red 
galaxies) since in this direction the BAO peak is better 
confined. 
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